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ABSTRACT 


A high-resolution, multi-level, primitive equation ocean model is used to examine 
the response of an idealized, flat-bottomed, eastern boundary oceanic regime on a beta- 
plane to both steady and daily-varying climatological wind forcing. The area of study 
is a coastal region within the Chile Current Svstem from 22°S to 34°S. When steady 
wind forcing is used, an equatorward surface current and poleward flowing undercurrent 
develop. Eddies are also generated, with initial formation in the poleward end of the 
domain. When daily-varying wind forcing is used, there is large spatial variability in the 
oceanic response. A relatively weak poleward flowing undercurrent appears, first in the 
poleward part of the domain. An equatorward surface current also develops and inten- 
sifies during the upwelling season. Eddies are generated and develop farther poleward 
in the domain than in the first experiment. The eddies are largest in the equatorward 
end of the domain. The eddy motion 1s closely tied to seasonal influences, with cyclonic 
(divergent) eddies traversing towards areas of higher dynamic heights when there are 
divergent wind fields present, and vice versa. In addition to the currents and eddies, 


upwelling and cold filaments are evident in both experiments. 
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I. INTRODUCTION 


The eastern boundary current system adjacent to the west coast of Chile, although 
well known for its biological productivity, has (until recently) attracted relatively little 
investigation into the physical oceanography of the region, despite an obvious relation- 
ship between the fertile waters in the region, and the fluid dynamics that make them so. 
The mean climatological current system is comprised of several major features, includ- 
ing: the equatorward flowing surface current called the Chile-Peru Current, and a 
poleward flowing countercurrent and/or undercurrent known as both the Gunther Cur- 
rent and the Peru-Chile Undercurrent. 

Eddies, jets and other instabilities that form along this and other eastern boundary 
current systems (e.g., the California Current System) have begun to generate a signif- 
icant amount of interest. Process-oriented studies using numerical models, such as those 
spearheaded by Batteen ev al. (1989), attempt to quantify the forces that create these 
features. 

Because the Chile Current Svstem has not been previously investigated with nu- 
merical models, the objective of this thesis will be to investigate the effects that wind 
forcing has on the generation of currents, eddies and filaments in the Chile Current sys- 
tem. A high-resolution, multi-level, primitive equation model will be used to simulate 
the effects of wind stress on the eastern boundary current region off Chile. The numer- 
ical model used for this study is one developed by Haney (1985), and adapted for eastern 
boundary current regions by Batteen (1989). 

The study is organized as follows. In Chapter II we review the oceanic and wind 
regime off the coast of Chile. In Chapter IJ] we discuss the numerical model used, as 
well as the specific initial and experimental conditions. In Chapter IV we describe the 
analysis techniques used to investigate the dynamical reasons for eddy generation. The 


modeling results are presented in Chapter V, and a summary 1s given in Chapter VI. 


II. BACKGROUND 


A. REGIONAL DESCRIPTION 
I. Surface current 

The mean climatological ocean adjacent to the coast of Chile is an eastern 
boundary current regime that consists of a slow-moving (~ 10 cm/s) surface current 
overlying a relatively strong (~ 25 cm/s) and steady coastal undercurrent (Wooster and 
Reid, 1963). The surface current is similar to other eastern boundary currents in that it 
is a sluggish, wide and generally shallow stream moving equatorward along a west coast 
(Wooster and Reid, 1963). Mass transports have been estimated to be between 
17.4 x 106 - 19.0 x 106 m/s (Wooster and Reid, 1963). The primary source of the Chile- 
Peru Current 1s, as expected, the West Wind Drift, because the South American conti- 
nent extends far enough south to intercept at least part of the flow, and divert it towards 
the north. Because the high latitudes are areas where precipitation typically exceeds 
evaporation, the northward flowing waters are relatively low in salinity, averaging less 
than 34.5 ppt. As the flow moves towards the north, salinity levels increase because of 
mixing and the general increase in evaporation along the way (Wooster and Reid, 1963). 

Tchernia (1980) and Gunther (1936) characterize the Chile Current System as 
consisting of two equatorward flowing currents, separated by a seasonal surface 
countercurrent. The outer current is usually referred to by those who adhere to this “two 
current’ philosophy as the Oceanic Chile-Peru Current; it extends to depths of never 
more than 1000 »: (Wooster and Reid, 1963), and usually no more than 700 m (Tchernia 
1980). The inner surface current is known as the Coastal Chile-Peru Current, with a 
maximum depth of ~ 200 m (Tchernia, 1980). 

Another description of the Chile-Peru Current 1s that it is a single, broad, slow, 
and sometimes intermittent surface current overlying a poleward flowing undercurrent 
(see Fonseca, 1989). In addition, Wooster and Reid (1963) describe the Chile-Peru 
Current as the principal component of the meridional circulation in the South Pacific. 

Yet another view of the Chile Current System is a combination of the previous 
discussions. Fonseca (1989) describes the current as being primarily one body of water; 
however, he states that in the southern domain the current 1s physically split into two 
components, due to either a countercurrent as described by Tchernia (1980), or to an 


undercurrent that has broached the surface (discussed in the next section). 


2. Peru-Chile Undercurrent 

Beneath the surface flow is a steady, poleward undercurrent, generally known 
as the Peru-Chile Undercurrent (PCU). The undercurrent 1s strong and steadv enough 
to be considered by some (e.g., Fonseca, 1989) to be the dominant ocean circulation 
feature in the Chilean eastern boundary current region, often broaching the surface of 
the ocean in spite of predominantly equatorward winds throughout the vear. Figure 1, 
from Silva and Fonseca (1983), illustrates the relationship between the surface current 
and undercurrent at different times of the year. The PCU 1s overlain by an equatorward 
surface current south of 30°S, but broaches the surface and widens north of 30°S. The 
PCU is wider and more of a dominant feature in the north during the summer months, 
and becomes narrower during the winter. In the south, the PCU is dominated by the 
surface current, especially in the summer months, when winds are more meridional in 
Orientation. Tchernia (1980) also acknowledges the existence of an undercurrent, de- 
scribing it as distinct and flowing beneath the countercurrent, with depths extending to 
several hundred meters. 

Until recently, the PCU had been relatively unexplored. Gunther (1936) dis- 
cussed the undercurrent that would eventually carry his name, but very little additional 
research was conducted on the PCU until 1961, when Wooster and Gilmartin (1961) 
measured it using parachute drogues. Their article sparked a new interest in the region, 
and soon several expeditions set out to study the waters adjacent to the South American 
coast. Between 1960 and 1982 the area was explored by twelve different expeditions 


(Fonseca, 1989). Results of the experiments revealed some consistent traits of the PCU: 
The PCU 1s present all year. 

It is easily distinguishable from the surface to about 600 m depth. 

Its core 1s at approximately 150-200 m depth. 

It is principally located between the coast and 200 km offshore. 


Average speeds are on the order of 25 cm/s. 


PD Oe eee 


Transport 1s estimated to be 10 x 10° 3/s. 


In a recent workshop on poleward flows along eastern ocean boundaries, 
Clarke (1989) concluded that, although many mechanisms may cause mean flows along 
eastern boundary regions, only two have thus far been shown to have a significant im- 
pact on the creation of a poleward flowing undercurrent. The first of these is the gen- 
eration of poleward flow by small-scale irregularities on the ocean bottom, in which 


coastally trapped lee waves are formed by the irregularities under the influence of a mean 


equatorward flow. Poleward flows of this type were estimated by Clarke (1989) to be 
of such a small magnitude (~ 2.5 cm/s) that we will neglect them in this study. The other 
mechanism mentioned by Clark (1989) that can generate poleward flows of significant 


magnitudes (~ 25 cm/s) 1s wind stress, the role of which will be studied 1n this paper. 


B. CLIMATOLOGICAL WINDS 

The Chile Current System is influenced predominantly by equatorward winds 
throughout the year. The direction, and to a large degree the magnitude of the winds 
are the result of a senu-permanent subtropical high pressure system that is similar in 
nature and behavior to its counterpart in the Northern Hemisphere, the North Pacific 
Subtropical High, as described in Nelson (1977). 

The South Pacific Subtropical High migrates meridionally with the seasons, reaching 
its northernmost extent during the austral winter, with its center located roughly at 
20°S. As the seasons progress, the high pressure system moves southward until its cen- 
ter reaches about 30°S during the summer months. 

Because of this migration, maximum values of wind stress and wind stress curl vary 
temporally at given locations. In Peru, for example, the maximum wind stress and wind 
stress curl values are seen during the austral winter months (Bakun and Nelson, 1991), 
presumably because the South Pacific Subtropical High reaches its northernmost mi- 
gration at that time, and the resulting winds take on more of an alongshore character- 
istic. The middle of the Chilean coastline, in contrast, does not experience maximum 
wind stress and wind stress curl values until late spring and summer, corresponding to 
the upwelling season in that region (Fonseca, 1989). Maximum wind stress curl values 
for the entire region are seen during the austral summer months, near 33°S. 

Winds stress curl values are usually at a relative minimum in the region off the 
northwest coast of Chile, compared to the rest of the Chilean coastal waters. Bakun and 
Nelson (1991) argue that this curl minima (i.e., anticyclonic curl) might be caused by the 
coastal winds blowing into the coastal bend near Arica, rather than remaining parallel 
to the coast. As the winds blow onto the land, alongshore speeds would in effect go to 
zero, eliminating any wind stress curl-induced upwelling. Wooster (1970) discusses a 
cyclonic cur] minimum in the same area, although he concentrates on the curl of the 
surface current, rather than on that of the wind. His discussion does, however, tend to 


Support the observation that cyclonic curl values are at a relative minimum near Anica. 





Figure 1. Surface current and undercurrent, from Silva and Fonseca (1983): Sim- 


plified representation of the Peru-Chile Undercurrent (shaded) between 
the coast of Chile and ~ 290 Am offshore. 


II. MODEL DESCRIPTION 


A. MODEL EQUATIONS 

The numerical model used in this study was originally developed by Haney (1985), 
and later adapted by Batteen ez al. (1989) for application to eastern boundary current 
regions. The following description of the model equations is based primarily on similar 
descriptions found in those papers. The multi-level, primitive equation (PE) model uses 
the hydrostatic, rigid lid, £-plane, and Boussinesq approximations. The governing 


equations are: 





Siren 
Equations of Motion: aE ra? =i = +f = A,V ut Kn = + 6<u) [1] 
zZ 
dye 
dt fo day fu 2] 
eu, 8 
u v 
v=-| ax Op \de [3] 
—H 
0 
0 0 
Hydrostatic Equation: 1) -| pede — — ({ pgde)dz [4] 
2 H _y Z 
Equation of State: Pp = po(I — a(T — T)) [5] 
72 
First Law: a =~ AT + Ky a ae OA!) [6] 


Important features of the above equations include: 


e A right-handed cartesian coordinate system is used in the above equations, with x 
pointing in the easterly direction, y pointing alongshore in a northerly direction, 
and z pointing upward. Accordingly, u 1s velocity in the x direction, and is positive 
when pointing toward the east; v 1S positive to the north, and w 1s positive upward. 


e T is temperature. 
°¢ pis density. 
e fis the Coriolis parameter. 


e¢ p’ is the departure from the vertically averaged pressure (which is assumed to be 
zero). Because the mean pressure is assumed to be zero, the barotropic mode ts 
considered to be ignored. 


® ¢€ 1S simply a constant of integration used in Equations (3) and (4). 


¢ Equation (5) assumes p to be a function of temperature only. Any dependence on 
salinity 1s ignored. 


e The terms 6{u), 6,(v) and 6{7) represent the vertical turbulent mixing of zonal and 
meridional momentum and heat, respectively, due to surface layer processes, by a 
dynamic adjustment mechanism. This adjustment maintains the dynamic stability 
in the water column. tn the sense of lighter water overlying denser water (Haney 
et al., 1978). 


The boundary conditions at the top (z = 0) of the model ocean are : 


Ama = 9 | BA 
Kn ge = 8] 
Ky Sh =- Qs (91 
w=0. [10] 


In Equation (8), t 1s the surface wind stress. In Equation (9), Q, is the net upward flux 
of longwave radiation. 


The bottom (z = —H) boundary conditions are as follow: 


moe = Crlu? + v*)(u cos y—vsin y) [11] 

K,. = Cp(u? + v?)'?(y cos y + usin y) [12] 
aT 

Spe [13] 

w=0 . [14] 


In the above equations, C, 1s a bottom drag coefficient and y = 10° 1s the geostrophic 
inflow angle (Weatherly, 1972). Table 1 provides a list of other symbols found in the 


model equations, as well as values of constants used throughout this study. 


B. MODEL DOMAIN AND RESOLUTION 
The domain of both experiments in this study consists of a rectangular area, with a 


meridional length of 1280 km, and a zonal width of 576 km. The eastern side of the box 


Table 1. VALUES OF CONSTANTS USED IN THE MODEL 
[[vaue =< Name SSOSCSCS~S~C~S 
Ta ada SS~S~*d Earth oration rate 
pC 0.958 calgm °K | Specific heat of sea water 
Te, x10" YY Atmospheric drag coefficient 
pL | S95 calgm' | Latent heat of sea water 
(Ti 278.2°K_ | Constant reference temperature 
pe | 23x10 gmoms | Density ofair 
— 1.0276 gm cnr Density of sea water at 7, 
po | 2.01 x 104 SK | Thermal expansion coefficient 
pK fT 0 Number of levels in vertical 
9.0 x 105 cm 
2.0 x 108 cm 
4.5 x 105 cn 
20x 10" ans 
2.010" ant 5 
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is assumed to touch the South American coast (Figure 2). In addition, for modeling 
purposes, the coastline 1s assumed to run exactly north-south; consequently, the domain 
is similarly oriented. The actual coastline is relatively straight, but runs in a slight 
southwest to northeast direction. At 34° S the western edge of the coast is approxi- 


mately 72° W; at 22° S, the coast is located at about 70° W. 


C. FINITE DIFFERENCE SCHEME 

The numerical model used in this study is a staggered grid model, using the Arakawa 
and Lamb (1977) B-scheme. Batteen and Han (1981) have shown that this scheme 1s 
appropriate when grid spacing 1s approximately on the same order as, or less than, the 
Rossby radius of deformation, which meets the criteria of this study. 

In the vertical, the model has ten layers, with constant z-levels of 13, 46, 98 182, 
316, 529, 870, 1416, 2283, and 3656 m. Consistent with Haney (1974), this spacing 
scheme concentrates more on the upper, dynamically active part of the ocean, above the 


thermocline. 


D. HEAT AND MOMENTUM DIFFUSION 

Because this study is designed to investigate the generation of mesoscale eddies 
along the coast of Chile, the tvpe of diffusion used in the model equations becomes an 
important issue. Laplacian lateral heat diffusion tends to suppress the formation of 
mesoscale eddies by diminishing baroclinic signals at that scale, thus making it more 
difficult for baroclinic instabilities to overcome diffusive damping (Holland and Batteen, 
1986). Biharmonic diffusion, on the other hand, tends to allow the formation of 
mesoscale instabilities, because it tends to damp at scales smaller than mesoscale eddies 
(Batteen er al., 1989). For these reasons, biharmonic lateral heat and momentum dif- 


fusion are used in this paper. Coefficients are listed in Table 1. 


FE. SURFACE THERMAL FORCING 

To isolate the effects of wind forcing in these process-oriented studies, the net heat 
flux at the sea surface should initially be zero. In this way, any heat flux that 1s gener- 
ated will have been generated by an increase or decrease of sea surface temperature re- 
sulting from wind forcing effects (Batteen et al., 1989). To accomplish this, surface 
generated heat flux (w'T’') is set equal to (S,) , which is defined as the solar radiation at 


the sea surface (Haney er a/., 1978). Simply stated: 


Sy — (w'T') = 0. eS] 


These terms and their components are in turn broken down into the following math- 


ematical definitions: 


Sp = So | PoC [16] 

(w'T "Jo = (Qp + Qs + Qe) | PoC [17] 

So = 0.95(0.74 — 0.6n)S; [18] 

S, = 3.642 x 107° [19] 

On = 0.9850(Ts) x [0.39 — 5.0(e,)'*](1 — 0.6’) [20] 
Os = paCuCy|VI(Ts — Ty) (21) 

Or = paCyL|V|(0.622/py4)(es — e4) [22] 


The variables used in these equations are defined in Table 2. Values for Ts, 7, py, ey, 
specific and relative humidity (used to compute e, ) are all taken from Esbensen and 
Kushnir (1981), while solar insolation values are taken from List (1984). To estimate 
annual means, the values were averaged over twelve months. Using these values, the 


ambient air temperature that results in zero heat flux at the surface 1s 290.13° K. 
F. BOUNDARY CONDITIONS 


The rectangular box surrounding the model domain 1s bordered on the northern, 
southern and western edges by open ocean. The eastern edge 1s assumed to border the 
western coast of Chile; it is a no-slip, straight vertical wall extending from the surface 


to 4500 m depth. At this eastern boundary, both tangential (v direction) and normal (x 


Table 2. VARIABLES USED IN COMPUTING NET HEAT FLUX 


[| value ——* «Description 


: Equation 19 Flux of solar radiation at top of 
atmosphere 


Ta, | Equation 20 | Net upward longwave rad 
Te | lebmejen? —SS~S~dCSWind ress SSS 
Co 0.958 calc! | Specheatofseawater 
C—O eS ae Fractional cloud cover 

PSH Reflectivity of sea surface 
ee x7 
1.757 x 107 
Te) 135x108 ~~~ Stefan Bottamann constant 
pe | 16.49 | Atmos vapor pressure 
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direction) components of velocity go to zero. A flat bottom at 4500 m is assumed 
throughout the entire domain. Both of the above assumptions can be justified for two 
reasons. First, because this 1s a process-oriented study, the inclusion of a rough 
coastline or bottom would induce additional shears and curl values (compared to those 
induced by a flat bottom and straight vertical coastline as in Batteen er al. (1989) and 


Fonseca (1989)), with the result that the particular process being investigated (i.e., wind 


forcing) would not be isolated. Second, the Chilean coastline is relatively straight and 
smooth, the continental shelf is relatively narrow, and the continental slope is steep, 
dropping from sea level to over 1000 fathoms (~ 2000 m) in a distance of less than 
25 km from the coast (Ziegler et al., 1957). Therefore, the use of a straight, vertical wall 
is NOt too unreasonable an approximation for a process-oriented study designed to iso- 


late the effects of wind forcing on the Chilean Current System. 


G. INITIAL CONDITIONS 

Climatological vertical and horizontal temperature profiles, derived from Levitus 
(1982), show that, over the study domain, the average seasonal surface temperatures 
vary between a summertime maximum of roughly 22° C at ~ 22°S, to a wintertime 
minimum of 13° C at ~ 34° S. The average annual sea surface temperature for the area 
is 17° C. Although a relatively weak north-south temperature gradient is evident at all 
depths, seasonal variations are negligible below 100 m depth. 

A temperature profile based on Levitus (1982) was used in both experiments to 
represent the initial temperature conditions for the model domain. The method for in- 
corporating the values of the profile into the model was similar to that used in Batteen 
and Rutherford (1990). Initially, mean annual temperatures were taken from Levitus 
(1982) for applicable 5° squares in or adjacent to the model domain, then interpolated 
in an east-west direction to match the middle of the domain. The temperature profile 
at each point was then interpolated vertically to correspond to the midpoints of the ten 
layers used in the model. These computations resulted in two “smoothed” temperature 
profiles, one representing the northern part of the domain (centered at 17.5° S), and the 
other representing the southern region (centered at 37.5°S ). These two profiles were 
interpolated once again, this time in a north-south direction, to obtain an annual tem- 


perature profile, centcred at 28° S. 
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As in Batteen er al. (1989, 1991) an exponential curve was then used to approximate 
the climatological temperature profile that had been computed from Levitus (1982). The 


form of the equation used was: 
T(z) = Ty + ATe7!". [23] 


The approximation uses a length scale of 450 m, and assumes 7, = 2° C to be the tem- 
perature at great depth. AJ = 15° C is the change of temperature between the bottom 
and top of the profile. Figure 3 shows that the resulting temperature profile fits the 
mean annual interpolated profile. This exponential profile was used to provide initial 


values at each of the ten levels in the model. 


H. EXPERIMENTAL CONDITIONS 
I. Experiment One: Constant Wind Stress 
In the first experiment, a constant equatorward wind stress of 1 dyne/cm? 
( ~ 830 cm/s) is applied to all gridpoints within the domain, except for a band at the 
northern and southern ends, where no wind forcing is applied. McCreary (1981) dem- 
onstrated that this technique allowed for the propagation of coastal Kelvin waves, which 
in turn allowed for the creation of an undercurrent and a coastal jet (Batteen er al., 
1989). Because the winds are predominantly equatorward all year off Chile, the model 
is run for 240 days, which is a sufficient time for currents and eddies to develop in the 
model. 
2. Experiment Two: Variable Climatological W inds 
In Experiment 2 , daily climatological winds are used to simulate the temporal 
and spatial wind variations that have been observed off the coast of Chile. Data, in the 
form of climatological wind stress fields, are taken from Bakun and Nelson (1991). The 
wind stress fields were derived from all available ship wind reports that have been col- 


lected in the region between 1850-1979, after eliminating physically unrealisti¢: reports. 


Figures 4 and 5, taken from Bakun and Nelson (1991), show the wind stress and wind 
stress curl fields, respectively, that are used. The wind stress fields are grouped into 
one-degree square regions (their midpoints being on the intersections between integer 
latitudes and Jongitudes), and a monthly mean value 1s computed for each region. 
Because wind speed values (vice wind stress values) are used in the model, wind 


speed values were obtained using the following equation: 


UF 7 








where 7 1s wind stress, p, is air density, C, is the atmospheric drag coefficent, and U is 
wind speed. The resulting wind values were then smoothed using a nonlinear 3 X 3 
median smoothing technique described in Rabiner er al. (1975), and used successfully in 
Batteen es al. (1991). This technique was shown in Rabiner er al. (1975) to preserve the 
Significant features in a field that can otherwise be hidden or “smoothed-out” when a 
linear filter 1s used. After smoothing, the resulting wind values for the 
12° latitude by 6° longitude region being studied were interpolated onto a 65 X 65 grid, 
basically covering the same geographical area. The easternmost gridpoints were con- 
sidered to be touching the continent. The resulting twelve monthly wind fields were then 
linearly interpolated into 365 daily fields, which were used to simulate one vear’s worth 
of daily, climatological wind values. In their final form, the data consisted of 365 grids 
(65 X 65), each gridpoint having a unique interpolated value. In the model, the appro- 
priate day’s values are read at the beginning of each integer day. An example of the 
wind fields can be seen in Figure 6, which shows the annual migration of maximum 
winds from the south in the austral summer, to the north in the austral winter months, 
the effects of which are discussed in Chapter V. 

It should be kept in mind that the daily wind values are a linear interpolation 


of monthly climatological winds, and not daily climatological values. As a result, the 


values should be considered as linear extensions of month-to-month trends, divided into 


dailv values to allow for a smooth transition from one monthly mean to the next. 
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Figure 2. Nlodel Domain: Area of study is a 576Am by 1280 km box off the 


coast of Chile. 


TEMPERATURE PROFILE 
ANNUAL 


Depth (m ) 





0 2 4 #6 8 10 122 4% 6 18 20 
Temperature (degrees Celsius ) 


Figure 3. Temperature profile used in both experiments: A smooth exponential 
curve (solid line) was used to approximate the climatological temper- 


ature profile (solid line with small circles) computed from 


Levitus (1982). 
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IV. ANALYSIS TECHNIQUES 


The following energy and spectral techniques, adapted from Batteen and Rutherford 
(1990), will be used in the following chapter to analyze the generation and stability of 


the Peru-Chile Current and eddies. 


A. ENERGY ANALYSIS 
An energy analysis based on that of Han (1975) and Semtner and Mintz (1977) 1s 
made to gain a better understanding of the energy transfers in the unstable flow. The 


energy calculations are presented using the Semtner and Mintz (1977) notation: 


() time average 

Ci time deviation 

() horizontal space average 
ei horizontal space deviation. 


The kinetic energy (K) 1s calculated by: 
a a [25] 


and presented in a time series plot. After reaching a quasi-steady state in which the total 


kinetic energy 1s nearly constant, the time mean and time eddy kinetic energy are calcu- 





lated by: 
—2) —2 
hie me [26] 
ary ame? 
=a = , [27] 


Zl 


Available potential energy (P) 1s calculated by: 


Pos] ry( 2) | (28) 


The temporal mean and eddy available potential energy are calculated by: 


P = | + Tr oe a | (29] 
P= | a. ¥( a yr (30] 


The transfers between the energy types are defined, after Semtner and Mintz (1977), 





(K + P} =— ag[Tw] [31] 
{P’ > K’} =agi/Tw’'] [32] 
Koxj=1(v.7v+ ivy) [33] 


(P= Py moe] Tv.wT"( SE) } [34] 


The model output consists of velocity components and temperature at daily intervals 
at each gridpoint and neither the vertical velocity nor any advection terms calculated by 
the model are stored. Because the calculation of the energy transfers requires both ver- 
tical velocity and numerous advection terms, these were recalculated in the same manner 


as in the model but using the stored values of u, v and T as the input data. Although 


Ze 


the energy transfers are recalculated from model output, they are consistent with the 
initial calculations of vertical velocity and advection terms in the model. 

Semtner and Mintz (1977) applied their energy transfer analysis to currents which 
had become unstable, generated eddies and then reached a quasi-steady state. In this 
study, the quasi-steady energetic state prior to eddy generation 1s examined, and the 
energy transfers are used to argue for the instability mechanism which leads to the sii 


eddy generation in each case. 


B. SPECTRAL ANALYSIS 

A spectral analysis technique is used to estimate the dominant wavelengths at which 
eddy growth occurs. The alongshore wavenumber ts selected based on an assumption 
of a meridional anisotropic preference for the eddy development. The model configura- 
tion of 65 alongshore temperature gridpoints and 64 velocity points is ideal in allowing 
data over the entire domain to be used to compute the discrete Fourier transform based 
on 64 points. Leakage due to the finite length of the series is reduced using a Hanning 


window. The loss of variance due to the windowing is then compensated by scaling the 


calculated one-sided spectral density appropriately. 


Ze 


V. RESULTS OF THE WIND FORCING EXPERIMENTS 


Experiment | studies the effects of steady, upwelling-favorable winds in the Chile 
eastern boundary current region. Experiment 2 investigates the effects of seasonally 
varying winds in the same region, using climatological monthly values that have been 


interpolated to dailv values. 


A. EXPERIMENT 1 

In Experiment 1, a constant wind stress was applied to the interior of the domain 
for 240 days. As in Batteen ef al. (1991), inertial oscillations of near-surface ocean cur- 
rents developed initially. After a few days of model simulation time, the oscillations were 
damped, leaving quasi-steady offshore Ekman transport, as shown by the surface veloc- 
ity vectors in Figure 7A. Figures 7B, 7C and 7D indicate that an equatorward-flowing, 
upwelling-favorable surface current had bcen generated. A poleward-flowing undercur- 
rent (not shown) had also developed. The surface current flowed at ~ 20 cm/s, and ex- 
tended ~ 50 km offshore, while the undercurrent, centered along the coast at ~ 300m 
depth, had weaker velocities of ~ 4 cm/s. 

The first evidence of developing eddies occurred by day 50, as evidenced by the ap- 
pearance of perturbations in the zonal velocity fields (Figure 8). (Perturbations are 
more clearly seen in these types of presentations than in meridional velocity plots, be- 
cause the latter tend to be dominated by the currents, which are predominantly in the 
meridional direction.) The time sequence of surface zonal velocity shows the develop- 
ment and growth of eddies from days 50-160. The baroclinic and barotropic transfers 
(based on Equations (32) and (33), respectively), are given in Figure 9, and show that 
both types of instabilities contribute to the development of eddies in the region. 


Barotropic instability is evident near the offshore region of the coastal jet and near the 
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coastline because of the horizontal shear in the region. Baroclinic instability, which 
draws its energy from the vertical shear in the mean flow (Haltiner and Williams, 1980; 
Batteen er al., 1991), was evident in the interior of the coastal region. 

By day 60, meanders were evident at y ~ 250 to 500 Am near the coast in the velocity, 
dynamic height and temperature fields, as seen in Figures 1OA-10D. As in the f-plane 
experiment of Batteen er al. (1989), eddies first formed in the poleward end of the region, 
as expected. Figures 10 and 1] show that the dominant wavelength at day 60 1s 
~ 100 km. 

Figure 12 shows three meridional velocity cross-sections averaged over days 50 - 60. 
The three cross-sections correspond to areas near the equatorward, middle, and 
poleward parts of the model domain, respectively. Figure 12 shows that the surface 
current has a width of ~ 100 Am, and an average velocity of ~ 18 cm/s, while the under- 
current has a core velocity of ~ 10 cm/s at ~ 250m depth. The undercurrent weakens 
towards the equator and intensifies near its poleward end. 

Figure 13 shows the fields of surface velocity vectors, temperature, velocity con- 
tours, and dynamic heights (relative to 2000 m) at day 120. The perturbations and 
meanders seen earlier in the southern part of the domain (see Figure 7) have developed 
into anticyclonic and cyclonic eddies. For example, a cvclonic eddv 1s discernible at 
y ~ 256 km with a corresponding cold core temperature of ~ 14.5°C. Figure 14 shows 
that, although both barotropic and baroclinic instability processes contribute to the 
formation and development of eddies, barotropic energy transfer 1s dominant, due to the 
larger energy transfer. Figure 15 indicates that the surface current has intensified to 

> 20 cm/s near the poleward end of the domain. The undercurrent maintains a core 
velocity of ~ 5 — 10 cm/s , and is better developed at ~ 620 km than towards the pole. 

Closed eddies were fully formed by day 160 (Figures 16A and 16D). While average 


values of the coastal surface current were ~ 30 cm/s (not shown), a maximum value of 
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~ 40 cm/s was reached by day 160 (Figure 16C). A cold, offshore-flowing filament ex- 
tends to ~ 350 km offshore (Figure 16B). Figure 17 shows that the current extends to 
~ 150 — 200 m depth in the equatorward end of the domain, while the undercurrent core 
is at ~ 350m and centered at ~ 10 km offshore (except at the poleward end, where the 
core is at the coast). Figure 18 shows that the dominant wavelength was ~ 250 km. 
The maximum velocity of the undercurrent was reached by day 210, when it reached 
a magnitude of ~ 25 cm/s (Figure 19). The core moves from ~ 10 km (Figure 19A) to 
~ 40 km (Figure 19C) offshore from north to south. What appeared to be a surface 
countercurrent in Figure 19 is most likely a manifestation of a cyclonic eddy that had 
formed in that area, with its poleward-flowing eastern side dominating the 
equatorward-flowing mean surface flow. This is supported by Figure 20D, which shows 
an upwelling-favorable low height center in the same region (y ~ 300 km), accompanied 
by a corresponding cold core center. Figures 20A and 20C both show poleward com- 
ponents in the vicinity of the eastern edge of the cyclonic eddy. Also evident are cold, 


offshore-flowing filaments extending as far as ~ 400 km offshore, as seen in Figure 20B. 


B. EXPERIMENT 2 

In Experiment 2, daily-varying wind forcing, corresponding to interpolated monthly 
climatological values taken from Bakun and Nelson (1991), were used in the model. The 
experiment ran from day 1 (corresponding to January |) to day 365. In the early part 
of the experiment, which corresponded to the middle of the austral summer, winds were 
generally meridional throughout the domain (see Figure 6), and the South Pacific 
anticyclonic gyre was at its southernmost extent. Figure 5, taken from Bakun and 
Nelson (1991), also shows a relative wind stress curl minima in the northeast corner of 
the domain during the summer months. This would imply that upwelling would be en- 
hanced in the southern portion of the domain, where cyclonic wind stress curl values 


were higher. As was the case in Experiment 1, an equatorward surface current was es- 
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tablished by day 20, extending ~ 80 km off the coast and reaching a magnitude of 

> 20 cm/s. An undercurrent first appeared at day 20 in the southern extreme of the 
region, but by dav 40 was present all along the coast, with a core ~ 10 km offshore, at 
~ 350 m depth (see Figure 21). 

A series of temperature cross-sections at day 50 (Figure 22) shows upward bending 
of isotherms near the coast, which is consistent with a more intense upwelling regime 
in the southern than in the northern part of the region. A planar view of the surface 
temperature at day 50 (Figure 23B) also shows a wider swath of cooler, upwelled water 
in the south, and a narrower band in the north. Surface velocity vector, zonal velocity, 
and dynamic heights in Figures 23A, 23C and 23D, respectively, show that, in contrast 
to Experiment 1, eddies had not yet formed by day 50, and velocities tended to be 
smaller, perhaps because of the weaker winds (generally less than 800 cm/s) throughout 
most of the region. 

The first evidence of perturbations was seen by dav 80 (Figure 24). Meanders occur 
~ 100 km to 200 km further to the north than their counterparts in Experiment 1. The 
reason for the difference between the two experiments was due to the different wind 
forcings used (i.e., uniform versus seasonally varying winds), because all other conditions 
remained the same. Spectral density analysis reveals a smaller wavenumber than had 
been present in Experiment 1, corresponding to a wavelength of ~ 85 Am (ieee Zo): 
(also seen in Figure 24, when eddies are being generated). Figure 26 shows that by day 
80 a weak, poleward flowing undercurrent 1s only evident in the middle and southern 
portions of the region, centered within 20 ky of the — at a depth of ~ 300m. Also 
shown 1s a surface current extending ~ 80 km off the coast, to a depth of ~ 150m. 

By day 140, the South Pacific High had migrated to the north, allowing a westerly, 
diffluent wind pattern to develop in the southern part of the domain (Figure 6), while 


winds in the north remained meridionally oriented. A time sequence of dynarhic height 


2] 


fields (Figure 27) shows that eddies were well developed by day 140, with larger eddies 
forming in the north, and smaller eddies in the south. As discussed in Batteen and 
Rutherford (1990), in a B-plane regime, barotropic instability tends to support eddies 
with radii greater than the Rossby radius of deformation, whereas baroclinic instability 
tends to support growth on the same order as the Rossby radius of deformation. As 
Figure 28 indicates, barotropic instability 1s prevalent through most of the coastal re- 
gions, while baroclinic instability tended to dominate only in the extreme south, below 
the region of diffluent winds. Because this effect was not seen in Experiment 1, the 
varying winds must have been the cause of the different sized eddies. By day 140 the 
swath of cooler, upwelled temperatures was no longer wider in the south (Figure 29). 
Instead, the widest band of cool temperatures was found in the north, once again likely 
due to the diffluent winds in the south, which tended to blow towards the coast rather 
than along it (Figure 6). In Figure 29, the time sequence of sea surface temperature 
contours also shows the development of cold filaments from days 80 to 140. 

By day 250, the winds had transitioned back into primarily a meridional regime 
(Figure 6). The time-averaged surface current had increased to > 50cm/s along the 
middle of the coast by day 250 (Figure 30), where wind magnitudes were greatest (Fig- 
ure 6). The dominant wavelength (Figure 31) had increased to ~ 250 4m, which cor- 
responded to the large eddies in the north. Smaller eddies were still evident in the south, 
and, as indicated in Figure 31, had become nearly as energetically significant as the 
eddies at longer wavelengths. The smaller scale eddy generation occurred during a pe- 
riod of primarily meridional, equatorward winds. 

Figure 32 shows a time sequence of surface dynamic height fields from days 250 to 
280, when a cyclonic, cold-core eddy at y ~ 700 Am breaks off and transits to the west. 
Figure 33 shows a time sequence of dynamic height fields towards the end of the year, 


when an anticyclonic, warm core eddy at y~ 768 Am travels eastward, into the region 
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of lower dynamic heights, possibly in response to a more meridionally oriented wind 
field. 

The equatorward surface current, which reached a maximum speed of ~ 50 cm/s by 
late austral spring (when winds had returned to a meridional! orientation throughout the 
region), remained relatively constant throughout the duration of the expermment. The 
average magnitude of the undercurrent (averaged over 30 days to filter out eddies) was 
~ 15 cm/s, at the end of the model simulation, when the South Pacific High had shifted 
back to the southern part of the domain. The largest instantaneous Ridcity of the 
undercurrent (not shown) was ~ 32 cm/s , occurring in the southern part of the region 
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Figure 7. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 50 of Experiment 1: Velocities are in cm/s, T is in °C, 


heights are in cm relative to 2000 m. Dashed lines are negative height 


values. 
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Figure 8. U (surface zonal velocity) for Experiment 1: Velocities are in cy/s. 


Dashed lines denote westward velocities. 
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Figure 9. Baroclinic and barotropic instability at day 50 of Experiment I. 
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Figure 10. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 60 of Experiment {: Velocities are in cen/s, T is in °C. 


heights are in cm relative to 2000 m. Dashed lines denote negative 
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Figure 11. Spectral density at day 60 of Experiment 1: Dominant wavelength is 


~ 100 km, where 1 = wavenumber" . 
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Figure 12. Zonal cross-sections of north-south velocity (V) at day 60 of Exper- 
iment 1: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 13. Surface current vectors, temperature, meridional velocity and dynatnic 
height at day 120 of Experiment 1: Velocities are in cm/s, T is in °C, 
heights are in cm relative to 2000 m. Dashed lines denote negative 


values. 
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Figure 14. Baroclinic and barotropic instability at day 120 of Experiment 1 
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Figure 15. Zonal cross-sections of north-south velocity (V) at day 120 of Exper- 
iment 1: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 16. Surface current vectors, temperature, meridional velocity and dynaniic 
height at day 160 of Experiment 1: Velocities are in cm/s, T is in °C, 
heights are in cm relative to 2000 m. Dashed lines are southward ve- 


locities (C) or negative values (D). 
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Figure 17. Zonal cross-sections of north-south velocity (V) at day 160 of Exper- 


iment I: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 18. Spectral density at day 160 of Experiment 1: Dominant wavelength 


is ~ 250 km. 
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Figure 19. Zonal cross-sections of north-south velocity (V) at day 160 of Exper- 
iment I: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 20. Surface current vectors, temperature, meridional velocity and dynamic 
height at day 210 of Experiment I: Velocities are in cm/s, T ts in °C, 
heights are in cm relative to 2000 m. Dashed lines are southward ve- 


locities (C) or negative values (D). 
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Figure 21. Zonal cross-sections of north-south velocity (V) at day 40 of Exper- 
iment 2: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 22. Zonal cross-sections of north-south temperature (T) at day 50: [em- 


peratures are in °C. 
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Figure 23. Surface current vectors, temperature, zonal velocity and dynamic height 
at day 50 of Experiment 2: Velocities are in cm/s, T is in °C, heights 
are in cm relative to 2000 m. Dashed lines are westward velocities (C) 


or negative values (D). 
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Figure 24. Surface zonal velocity (U) for Experiment 2: Velocities are in cm/s. 


Dashed lines are westward velocities. 
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Figure 25. Spectral Density at day 80 of Experiment 2: Dominant wavelength 


is ~ 250 km, with a secondary wavelength evident at ~ 85 km. 
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Figure 26. Zonal cross-sections of north-south velocity (V) at day 80 of Exper- 
iment 2: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 27. Dynamic height, days 80-140 of Experiment 2: Heights are in cm rel- 


ative to 2000 m. Dashed lines denote negative height values. 
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Figure 28. Baroclinic and barotropic instability at day 140 of Experiment 2. 
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Figure 29. Surface temperature contours for Experiment 2: Temperatures are 
in °C. Temperature bands in the south are narrower during the austral 


fall. 
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Figure 30. Zonal cross-sections of north-south velocity (V) at day 250 of Exper- 
iment 2: Velocities are in cm/s. Southward flowing currents are 


dashed. 
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Figure 31. Spectral density at day 250 of Experiment 2: Dominant wavelength 


is ~ 250 km, with a secondary wavelength evident at ~ 85 kon. 
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Figure 32. 
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VI. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

This study used a high-resolution, primitive equation numerical model to study the 
effects of steady and daily varying winds on the Chile Current System. The model do- 
main was a rectangular region extending 6° off the coast, between 22° and 34° Se. 

In Experiment 1, a constant wind stress was applied to the model domain, excluding 
the northern and southern open boundaries, for a period of 240 days. The objective was 
to simulate the generally equatorward-flowing winds that are seen in that part of the 
world, for most of the year. Both an equatorward-flowing surface current and a 
poleward-flowing undercurrent developed. Average surface equatorward and subsurface 
poleward velocities were ~ 25 — 30 cm/s. Eddies formed in the southern part of the do- 
main first, then spread north and south, eventually dominating the coastal region. The 
dominant wavelength for the eddies was ~ 250 km. 

In Experiment 2, temporally and spatially interpolated values of monthly 
climatological winds were used as the forcing for the model. Unlike Experiment 1, the 
varying winds allowed the following processes to be studied. First, it allowed the study 
of the effects of seasonal wind variations on the behavior of the current system. Second, 
it allowed the influence of wind stress curl on coastal dynamics to be investigated. Third, 
it introduced the zonal component of the wind stress, allowing the study of the oceanic 
response to this additional component of stress. 

Initially, when winds were generally equatorward, the two models were similar, with 
the development of coastal currents and upwelling along the coast. As the seasonal re- 
gime of the winds changed, differences, as expected, were discernible. For example, the 


poleward-flowing undercurrent in Experiment 2 was slower to develop, especially in the 
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northern parts of the region. The surface equatorward current reached a maximum 
speed of > 50 cm/s in Experiment 2 compared to a maximum speed of ~ 30 cm/s in Ex- 
periment 1. The average undercurrent in Experiment 2 was ~ 15 cm/s, which is slower 
than in Experiment | (which had velocities of ~ 25 cm/s ). In Experiment 2, the under- 
current was centered at ~ 300 m depth, compared to ~ 250 m in Experiment 1. In Ex- 
periment 2, instabilities did not appear until day 80, a full 30 days later than in 
Experiment 1. Baroclinic instability played a smaller role (than in Experiment 1) in the 
creation of eddies, except in the extreme southern part of the region. Once eddies did 
form, however, their sizes varied with latitude and time, tending to grow larger in the 
north where the winds were usually mostly meridional, and virtually disappearing in the 
south, when the alongshore wind stress was at a minimum. These southern eddies began 
to grow again, however, once stronger winds returned to the region. In the north, the 


dominant wavelength was ~ 250 km, while in the south it was ~ 85 km. 


B. RECOMMENDATIONS 

Of critical importance to anv process-oriented study is the comparison of model 
output with available observations. Without such comparisons, model studies are diffi- 
cult to verifv, and become speculative in nature. As mentioned bv several authors (e.g., 
Fonseca, 1989; Bakun and Nelson, 1991), the waters off the coast of Chile are 
notoriously lacking in both wind and current data, and until more observational studies 
are conducted, the Chile Current System will remain somewhat of a mystery. Johnson 
(1989) suggests that satellite imagery of the South American coastal waters may provide 
the most readily obtainable basis for increasing the number of available observations. 
Much of the model output from numerical simulations could then be compared directly 
to the imagery. Figure 34A shows a satellite image from Fonseca (1989) of a cold fila- 
ment extending out ~ 40 km off the coast of Chile, near 33° S, while Figure 34B shows 


the surface temperature field at day 290 of the Experiment 2 model run, in which cold 
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filaments can be seen extending to the west. The 14° contour at several locations 
protrudes 40-60 km from the coast, a distance very similar to that seen in satellite im- 
agery. Comparison of the two figures suggests that the model results may have some 
validity, although many more comparisons would need to be made before such a claim 
could be made with anv certainty. This merging of model output and real-time satellite 
verification may well be the next logical step in understanding the Chilean coastal wa- 
ters, and is a recommended area of future study, as such observations (e. g., GEOSAT, 
SEASAT) become available. 

The knowledge gained in this investigation can benefit Naval Oceanography bv 
providing a basis for prediction not only in the Chilean Current System, but also in other 
dynamically similar Eastern Boundary Current regions. An understanding of the for- 
mation and evolution of eddy and jet structures is an essential step in the development 
of a real-time analysis and prediction capability for these important features (Robinson 


et al., 1984). 
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Temperature contours from (A) _ satellite and (B) model 
output: (A) Satellite-derived temperatures from f[fonseca (1989); 
(B) Temperature contours from day 290 of Experiment 2. All tem- 
peratures are in °C. The inside box at y ~ 256 km in (B) corresponds 


to the size of the region shown tn (A). 
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